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annular  effect.  To  the  best  of  the  authors'  knowledge  these  are  the  first  solutions  demonstrating  either 
Richardson's  annular  effect  or  refraction  effects  that  are  obtained  through  solution  of  the  complete  Navier- 
Stokes  equations. 

To  help  discriminate  between  acoustic  energy  loss  due  to  artificial  numerical  dissipation  and  energy  loss 
due  to  physical  processes  (including  acoustic  energy  transfer  to  the  mean  flow)  a  second  study  was  conducted  in 
which  acoustic  wave  propagation  in  a  tube  with  no  mean  flow  was  investigated.  The  results,  in  agreement  with 
experimental  data,  demonstrated  that  mean  flow  is  excited  by  the  acoustic  flow  field,  a  phenomenon  referred  to 
as  acoustic  streaming.  Acoustic  energy  dissipation  was  significantly  lower  (per  wave  length)  for  this  test  case 
than  for  the  acoustic  refraction  test  case,  indicating  that  although  artificial  energy  dissipation  due  to  numerical 
errors  may  not  be  completely  eliminated,  this  energy  loss  is  small  in  comparison  to  acoustic  energy  dissipation 
due  to  physical  processes. 

The  third  study  investigated  acoustic  wave  propagation  in  a  tube  with  a  coexisting  mean  flow,  where  the 
mean  flow  is  injected  into  the  tube  through  its  lateral  boundary.  Results  demonstrate  the  existence  of 
Richardson's  annular  effect,  the  growth  of  the  acoustic  boundary  layer  due  to  wall  blowing,  the  spatial  growth 
of  the  pressure  and  axial  velocity  mean  shifts  (i.e.  growth  of  energy  contained  in  the  mean  flow  field  by  energy 
transfer  from  the  acoustic  flow  field),  and  the  formation  of  a  large  new  downstream  directed  axial  velocity  zone 
at  the  edge  of  the  acoustic  boundary  layer,  a  phenomenon  never  reported  before.  Acoustic  energy  losses  and 
mean  pressure  and  axial  velocity  shifts  were  significantly  higher  for  this  test  case  than  for  the  previous  two  tests 
due  to  the  addition  of  flow  turning  losses.  Though  the  present  research  was  able  to  verify  the  existence  of  flow 
turning  loss  phenomenon,  more  theoretical  and  experimental  research  efforts  are  needed  to  help  quantify  flow 
turning  loss  and  the  several  flow  phenomena  observed  in  this  investigation. 
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NUMERICAL  STUDY  OF  FLOW  TURNING  PHENOMENON 


INTRODUCTION 


The  growth  or  decay  of  random  small  amplitude  pressure  oscillations  in  combustion 
systems  is  determined  by  a  delicate  balance  between  the  driving  and  damping  (i.e.,  sources 
and  sinks)  of  oscillatory  energy  in  the  combustor.  When  the  driving  mechanisms  outweigh 
the  damping  effects,  small  amplitude  oscillations  that  are  indigenous  to  the  combustor 
chamber  (i.e.  random  noise  produced  by  the  combustion  or  turbulence)  can  be  amplified;  a 
condition  known  as  combustion  instability.  In  solid  or  liquid  propellant  rocket  motors,  these 
sustained  oscillations  may  affect  guidance  systems  and  alter  the  thrust  time  characteristics 
of  the  motor,  or,  in  extreme  cases,  may  result  in  significant  overpressures  and  motor  case 
failure. 

The  determination  (both  by  theoretical  predictions  and  experimental  measurements) 
of  the  various  acoustic  energy  sources  and  sinks  in  the  combustion  chamber  is  critical 
to  our  ability  to  predict  a  priori  the  stability  characteristics  of  proposed  solid  propellant 
motor  designs.  The  processes  which  dissipate  acoustic  energy  include  convection  and  ra¬ 
diation  of  acoustic  energy  through  the  nozzle  (nozzle  damping),  the  viscous  and  thermal 
losses  produced  by  interaction  of  the  condensed  phase  combustion  products  with  the  com¬ 
bustion  gases  (particle  damping),  the  nonlinear  viscoelastic  characteristics  of  the  grain 
and  case  (structural  damping),  viscous  losses  in  the  gas  phase  and  at  inert  surfaces,  gas 
phase  vibrational  relaxation  effects,  radiation  of  energy  through  the  motor  case  and  in¬ 
elastic  acceleration  of  combustion  products  leaving  the  propellant  surface  (flow  turning). 
The  primary  source  of  oscillatory  energy  considered  to  date  has  been  the  response  of  the 
propellant  combustion  zone  to  acoustic  pressure  and  acoustic  velocity  oscillations  (termed 
pressure  and  velocity  coupling,  respectively).  Other  sources  of  energy  include  the  com¬ 
bustion  of  incompletely  reacted  products  in  the  chamber,  distributed  combustion  of  metal 
droplets  in  the  gas  phase,  and  conversion  of  mean  flow  energy  to  acoustic  energy  due  to 
vortex  shedding.  The  objective  of  the  present  research  project  is  to  seek  an  understanding 
of  the  physical  mechanisms  by  which  energy  is  exchanged  between  the  mean  and  acoustic 
flow  fields  in  resonant  combustion  chambers  (in  particular,  solid  rocket  motors),  processes 
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that  have  been  the  subject  of  much  speculation  in  the  past,  but  which  have  never  been 
thoroughly  investigated.  Specifically,  the  phenomena  to  be  investigated  are  the  interaction 
between  the  acoustic  field  in  a  chamber  and  a  free  shear  layer  (i.e.,  vortex  shedding)  and 
the  so  called  “flow  turning”  effect. 

FLOW  TURNING 

The  term  flow  turning  is  commonly  used  to  describe  the  loss  of  acoustic  energy  by  the 
acoustic  field  in  a  combustor  resulting  from  the  inflow  of  combustion  products  through  the 
lateral  boundary  of  a  combustion  chamber  containing  longitudinal  acoustic  waves.  The 
mean  flow  of  solid  propellant  combustion  products  entering  normal  to  the  surface  must 
turn  to  the  axial  direction  parallel  to  the  boundary  in  order  to  exit  through  the  nozzle; 
simultaneously,  each  fluid  element  of  the  mean  flow  must  acquire  acoustic  energy.  It  is 
this  second  process,  involving  the  irreversible  action  of  the  unsteady  field  in  the  chamber 
upon  the  entering  flow,  which  has  been  hypothesized  to  be  a  loss  of  energy  for  the  existing 
acoustic  field.  This  energy  loss  is  often  referred  to  as  flow  turning  loss  because  mean  flow 
turning  is  necessarily  involved.  However,  it  is  really  a  rate  of  energy  loss,  proportional  to 
the  velocity  of  the  incoming  flow  perpendicular  to  the  boundary. 

The  modeling  of  flow  turning  involves  multi-dimensional,'  rotational  and  viscous  flow 
effects  on  widely  varying  distance  scales.  These  effects  only  increase  the  difficulty  of  model¬ 
ing  this  problem  due  to  difficulties  in  both  analytical  and  numerical  procedures.  Lack  of  a 
basic  fundamental  theory  of  mean  flow-acoustic  interactions  has  resulted  in  disagreement 
among  researchers  about  a  very  basic  point;  namely,  is  flow  turning  a  surface  effect,  a 
volume  effect  or  a  combination  of  the  two?  The  experimental  measurement  of  this  acous¬ 
tic  energy  loss  mechanism  is  difficult  since  this  effect  usually  appears  in  conjunction  with 
other  gain  and  loss  mechanisms  and  cannot  be  easily  isolated.  The  pioneering  theoretical 
work  of  Culick1'2  utilizing  a  one-dimensional,  inviscid  approximation  is  qualitative  at  best. 
The  logic  behind  these  approximations  was  challenged  by  several  researchers  working  in 
the  field  3’4.  Flandro’s  model3  is  locally  two-dimensional  and  laminar.  However,  Flandro 
chose  to  neglect  effects  of  rotational  core  flow  and  mean  flow  viscosity  and  to  limit  the 
model  to  small  Mach  numbers  and  small  disturbance  amplitudes.  Chung’s5  formulation 
utilizes  a  three-  dimensional,  viscous  analysis.  The  result  of  this  analytical  approach  has 
not  yet  been  compared  to  experiments.  Culick6  has  qualitatively  demonstrated  the  exis¬ 
tence  of  this  flow  turning  phenomenon,  but  due  to  the  very  low  Mach  number  facility  used 
to  conduct  the  experiment,  quantitative  information  could  not  be  obtained. 


Hersh7  conducted  an  experimental  investigation  of  the  flow  turning  losses  associated 
with  the  injection  of  steady  lateral  flow  of  cold  air  into  a  rectangular  duct  containing 
longitudinal  acoustic  waves.  This  experiment  was  designed  to  simulate  injection  of  hot 
combustion  products  into  a  combustion  chamber.  The  results  indicated  that  the  flow  turn¬ 
ing  process  absorbs  sound  (i.e.,  damping  of  acoustic  energy).  A  one-dimensional  model 
developed  by  Hersh8  seriously  underpredicted  the  measured  acoustic  energy  losses.  Hersh 
speculated  that  in  addition  to  the  flow  turning  absorption,  there  is  another  energy  ab¬ 
sorption  mechanism  related  to  acoustic  refraction.  Due  to  refraction  by  the  mean  flow 
velocity  gradients,  the  sound  pressure  near  the  side  wall  increases,  thereby  enhancing  the 
acoustic  energy  absorption  by  the  finite  admittance  walls.  Nevertheless,  Hersh  concluded 
that  "most  of  the  sound  absorption  appeared  to  take  place  within  the  fluid  interior  in 
support  of  Culick’s  ideas,  but  final  interpretation  of  the  data  requires  the  development  of 
a  two-dimensional  model  of  the  flow  turning  process". 

It  is  important  to  recognize  that  in  addition  to  flow  turning  and  refraction,  several 
other  physical  mechanisms  may  have  to  be  considered  when  an  attempt  is  made  to  properly 
model  acoustic-mean  flow  interactions.  It  had  been  demonstrated  through  analysis9-12  and 
experiments13-10  that  the  interaction  between  the  mean  and  acoustic  fields  may  (a)  mod¬ 
ify  the  turbulent  structure  of  the  flow;  (b)  accelerate  transition  to  turbulence;  (c)  cause 
flow  reversal  during  a  portion  of  the  wave  cycle;  and  (d)  inhibit  or  augment  heat  trans¬ 
fer  to  the  boundaries.  It  has  been  demonstrated1'  that  oscillating  flows  in  a  tube,  even 
without  mean  flow,  exhibit  fully  turbulent  behaviour  for  Reynolds  numbers  higher  than 
3000.  Richardson918  measured  the  root  mean  square  (termed  here  rms)  velocity  field  in  a 
circular  tube  in  the  absence  of  mean  flow  and  identified  the  overshoot  of  the  axial  veloc¬ 
ity  at  the  edge  of  the  acoustic  boundary  layer,  a  phenomenon  now  commonly  referred  to 
as  Richardson's  annular  Effect.  Mohajery19  and  BogdanofF1'*  observed,  utilizing  hot-wire 
measurements,  that  eddy  viscosity  increases  significantly  near  the  wall  for  acoustic  wave 
propagation  in  the  presence  of  mean  flow  at  Reynolds  numbers  up  to  10r>.  More  recently, 
Clamen15  conducted  a  series  of  experiments  using  an  hydrogen  bubble  technique  to  study 
water  flow  in  a  vibrating  tube  for  low  Reynolds  numbers  for  a  range  of  amplitudes  and  fre¬ 
quencies.  Among  the  theoretical  analyses.  Sexl20  obtained  a  solution  for  wave  propagation 
in  a  tube  with  no  mean  flow,  Ucliida10  extended  his  results  to  include  laminar  mean  flow, 
and  Romie11  investigated  harmonic  wave  propagation  in  laminar  pipe  flows  for  the  study 
of  heat  transfer  mechanism.  More  recently.  Barnett-’1  investigated  the  effect  of  longitudi¬ 
nal  pressure  oscillations  n  the  velocity  distribution  in  laminar  and  fully  developed  pipe 
flows.  Barnett’s  analysis  solved  the  incompressible,  Reynolds-averaged,  linearized  Navier- 
Stokes  equations.  The  solution  was  applied  only  to  low  Reynolds  number  flows.  Barnett 


demonstrated  the  existence  of  Richardson'  annular  effect  and  showed  that:  (a)  acoustic 
boundary  layer  thickness  is  strongly  influenced  by  mean  flow  Reynolds  number  and  by 
the  frequency  of  oscillations.  For  frequencies  of  interest  in  this  research,  but  for  signifi¬ 
cantly  lower  Reynolds  numbers,  the  analysis  predicted  (based  on  a  model  developed  by 
Richardson22)  the  acoustic  boundary  layer  thickness  to  be  within  0.00G5  radii  (for  1000  Hz) 
to  0.002G  radii  (for  G000  Hz);  (1))  acoustic  boundary  layer  thickness  and  maximum  axial 
velocity  overshoot  vary  during  a  cycle;  and  (c)  the  overshoot  at  the  edge  of  the  acoustic 
boundary  layer  varies  directly  with  frequency. 

To  properly  resolve  energy  transfer  between  the  acoustic  and  mean  flow  field,  it  is 
necessary  to  utilize  a  numerical  integration  scheme  that  has  very  small  diffusive  and  dis¬ 
persive  errors.  In  addition,  the  numerical  scheme  chosen  has  to  handle  very  fine  grid 
resolution  near  the  wall  in  order  to  resolve  energy  transfer  within  the  acoustic  boundary 
layer,  and  a  coarse  grid  resolution  near  the  centerline.  These  grid  variation  requirements 
can  be  best  satisfied  by  utilizing  implicit  integration  schemes2'1.  A  drawback  of  the  implicit 
codes  is  that  they  often  utilize  artificial  viscosity  to  stabilize  the  solution.  In  addition,  the 
truncation  and  splitting  errors  of  the  algorithms  may  act  as  artifical  energy  dissipation 
mechanisms.  These  numerically  induced  energy  dissipation  mechanisms  can  be  minimized 
by  reducing  the  artificial  viscosity  coefficient  value  used  in  the  calculations  and  reducing 
the  numerical  integration  time  step.  Nevertheless,  these  erroneous  dissipation  mechanisms 
cannot  be  eliminated  entirely. 

The  primary  objective  of  this  research  is  to  calculate  the  so  called  flow  turning  loss. 
This  is  done  by  calculating  the  axial  and  radial  acoustic  energy  distributions  for  a  down¬ 
stream  propagating  wave  as  it  traverses  the  length  of  a  tube  in  which  the  mean  flow 
is  introduced  into  the  tube  through  its  perimeter.  Before  calculating  acoustic  energy  it 
is  necessary  to  distinguish  between  acoustical  properties  (properties  that  travel  with  the 
speed  of  sound)  from  properties  that  are  convccted  by  the  flow.  The  separation  of  these 
properties  in  a  flow  that  involves  multi-dimensional  rotational  and  viscous  flow  effects  has 
never  before  been  extensively  studied  or  resolved.  The  most  important  question  namely, 
the  exact  definition  of  acoustic  energy  in  a  rotational  viscous  flow,  is  left  unanswered  for 
the  time  being.  Previous  analytical  work  on  this  subject  is  currently  beeing  evaluated 
and  its  application  to  this  specific  problem  will  be  made  in  the  future.  In  the  meantime, 
acoustic  properties  will  be  defined  here  as  the  difference  between  the  time  dependent  and 
the  steady  properties.  In  order  to  get  an  approximate  idea  about  acoustic  energy  losses 
in  the  system,  the  concept  of  acoustic  energy  that  was  derived  for  plane  wave  propagation 
in  irrotational,  inviseid  flow  (i.e.  the  time  averaging  of  pw)  was  adopted.  This  paper  will 
concentrate  on  the  description  of  the  computational  results  obtained  to  date.  Analysis  of 


4 


the  separation  of  the  calculated  properties  and  exact  calculation  of  acoustic  energy  transfer 
will  be  the  subject  of  future  research. 

The  first  case  to  be  described  is  acoustic  wave  propagation  in  a  tube  with  a  coexisting 
sheared  mean  flow.  Acoustic  energy  was  expected  to  diminish  due  to  both  acoustic  energy 
transfer  to  the  mean  flow,  and  through  artificial  energy  dissipation  due  to  the  numerical 
errors  of  the  integration  scheme.  To  help  quantify  acoustic  energy  dissipation  due  to 
numerical  errors,  the  calculations  were  repeated  with  significantly  more  mesh  points  both 
in  the  radial  and  axial  directions  simultaneously  with  a  reduced  integration  time  step. 
Acoustic  wave  propagation  in  a  tube  without  mean  flow  will  be  described  next.  This 
test  case  was  investigated  to  provide  a  baseline  for  evaluating  the  amount  of  erroneous 
numerical  dissipation  present  in  the  solution.  Though  Richardson's  annular  effect  and 
acoustic  streaming  had  been  observed  for  acoustic  wave  propagation  with  no  mean  flow18, 
it  was  expected  that  for  this  test  case,  energy  u ansfor  between  the  acoustic  and  mean  flows 
would  be  significantly  reduced  and  most  of  the  acoustic  energy  dissipated  would  result  due 
to  numerical  errors. 

The  third  test  case  to  be  described  will  be  acoustic  wave  propagation  in  a  tube  with  a 
coexisting  sheared  mean  flow,  where  the  mean  flow  was  injected  into  the  chamber  through 
part  of  its  lateral  boundary.  Comparison  of  axial  and  radial  acoustic  energy  distributions 
for  the  three  cases  investigated  should  enable  us  the  determination  and  quantification  of 
the  various  acoustic  energy  dissipation  mechanisms. 


ACOUSTIC  REFRACTION 


^  lion  sound  waves  propagate  in  a  moving  medium  (i.e.,  with  a  coexisting  mean  flow), 
the  main  physical  processes  affecting  its  propagation,  in  the  absence  of  viscosity  and  ther¬ 
mal  conduction,  are  convection  and  refraction.  The  sound  is  convected  by  the  moving 
medium  at  its  local  velocity.  In  the  presence  of  a  sheared  mean  flow,  local  convection 
velocities  are  a  function  of  position.  This  results  in  pressure  gradients  in  the  direction 
normal  to  the  flow  direction  (i.e..  in  the  radial  direction,  for  the  case  of  flow  propagation 
in  a  hard-wall  tube),  which,  in  turn,  results  in  the  excitation  of  velocities  normal  to  the 
flow.  Thus,  for  the  case  of  flow  in  a  rigid  tube,  when  the  sound  wave  travels  with  the 
mean  flow  (downstream),  the  effect  of  velocity  gradients  is  to  refract  the  sound  toward  the 
wall.  When  the  sound  propagates  against  the  mean  flow  (upstream)  the  effect  of  radial 
velocity  gradient  is  to  refract  the  sound  toward  the  centerline.  Schematics  of  downstream 
and  upstream  sound  propagation  are  shown  in  Figs,  la  and  lb,  respectively. 

Previous  research  efforts  to  model  and  understand  acoustic  refraction  phenomenon 
utilized  linearized  models2'1-2 ‘ .  For  initially  planar  acoustic  waves  propagating  down¬ 
stream  in  a  hard  wall  duct,  the  calculation  indicated  that  the  acoustic  pressure  at  the  wall 
is  100  db  larger  than  at  the  centerline.  This  study  also  indicated  that  acoustic  refraction 
effects  increase  with  mean  flow  velocity  and  with  the  frequency  of  the  oscillations.  In 
succeeding  studies20-2'  ,  the  linearized  treatment  has  been  extended  to  treat  sound  prop¬ 
agation  upstream  in  a  long  rigid  tube  as  well  as  the  effect  of  varying  shear  layer  thickness. 
These  studies  indicated  that  for  centerline  Mach  numbers  in  the  range  of  0.05  to  0.10  (flow 
velocities  which  are  typical  in  solid  motor  chambers),  the  acoustic  pressure  amplitude  near 
the  propellant  surface  may  be  twice  as  large  as  the  pressure  amplitude  at  the  centerline. 
Although  nonlinear  effects  can  be  expected  to  significantly  reduce  the  actual  refractive  en¬ 
hancement  of  the  near  wall  pressure,  these  results  indicate  that  acoustic  refraction,  which 
has  been  virtually  ignored  in  solid  rocket  stability  analyses,  may  actually  have  quite  im¬ 
portant  effects.  This  results  from  the  fact  that  both  the  steady  and  transient  response  of 
the  burning  solid  propellant  are  dependent  upon  the  amplitude  of  the  pressure  and  velocity 
oscillations  near  the  burning  surface. 

An  experiment  to  measure  acoustic  refraction  (-fleets  is  currently  in  progress8,28.  In 
this  experiment  air  is  blown  through  a  long  rigid  tube  in  which  a  plane  wave  is  excited  at 
the  inlet.  The  mass  flux  through  the  tube  was  varied,  resulting  in  a  centerline  Mach  number 
variation  from  0.03S  to  0.10.  The  results,  to  date,  indicate  that  the  acoustic  pressure  is 
amplified  near  the  pipe  wall:  the  degree  of  amplification  increases  with  Mach  number. 


With  a  centerline  Mach  number  of  0.1  (Re  =  2.7  x  lOM  the  measurements  indicated  a 
near  wall  acoustic  pressure  amplification  of  40%.  Such  an  increase  in  acoustic  pressure 
amplitude  may  be  critical  for  motor  stability  analysis.  The  discrepancy  between  the  linear 
theory  and  the  experimental  data  may  possibly  result  from  nonlinear  mean  flow-acoustic 
field  interactions,  which  will  tend  to  damp  acoustic  pressure  growth  near  the  walls.  These 
processes  arc  not  modeled  by  the  linear  analysis. 

NUMERICAL  SOLUTION 

Under  the  current  research  program,  the  MINT  code  was  acquired  and  modified  for 
the  investigation  of  acoustic  refraction  and  flow  turning.  The  code  solves  the  compressible 
Reynolds-averaged  Navier-Stokes  equations  using  an  efficient,  non-iterative  time  depen¬ 
dent  Linearized  Block  Implicit  (LBI)  scheme.29-30.  The  current  model  utilizes  either  a 
turbulent  kinetic  energy-algebraic  length  scale  (k  —  /)  model31  or  a  transport  turbulence 
(K  —  e)  model32.  The  later  includes  modifications33  that  allow  the  K  —  e  model  to  be  utilized 
throughout  the  entire  viscous  sublayer  without  any  wall  function  assumptions.  The  result¬ 
ing  code  was  verified  by  comparing  the  predictions  to  experimental  data34  for  flat  plate 
boundary  layer  flows  with  and  without  wall  transpiration.  For  a  complete  description  of 
the  equations  solved  and  the  turbulence  models  used  see  Reference  33. 

The  study  of  acoustic  refraction  utilizes  a  tube  with  a  radius  of  0.05  m  (two  inches) 
and  a  length  of  approximately  five  acoustic  wavelengths  (for  any  frequency  studied).  The 
temperature  and  the  axial  and  radial  components  of  the  velocity  are  specified  at  the  inlet 
plane,  while  constant  static  pressure  was  specified  at  the  outflow  plane.  The  calculation  is 
initiated  with  a  one-seventh  power  law  velocity  profile  at  the  inlet  plane,  having  a  centerline 
Mach  number  of  0.1.  After  convergence  to  steady  state,  the  pressure  and  axial  velocities 
at  the  upstream  boundary  are  perturbed  continously  to  produce  a  system  of  downstream 
traveling  waves  in  the  computational  domain.  An  exhaustive  study  of  non- reflective  bound¬ 
ary  conditions  demonstrated  the  inadequacy  of  all  currently  available  models  for  low  Mach 
number  sheared  flows.  Thus,  the  solution  was  terminated  when  the  front  of  the  initial 
wave  started  reflecting  from  the  downstream  boundary. 

Initial  results  for  the  refraction  study  were  presented  in  Reference  23.  These  results 
were  calculated  on  the  Air  Force  Weapons  Laboratory  Cray  1-S  computer.  Many  of  these 
calculations  were  repeated  on  the  Naval  Research  Laboratory  Cray  X-MP  computer  that 
has  twice  the  online  memory  storage  (core)  and  twice  the  speed  of  the  Cray  1-S.  The  repeat 
calculations  were  conducted  with  approximately  twice  as  many  mesh  points  per  wave 
length;  a  larger  total  number  of  points  in  the  axial  direction  (currently  about  250)  which 
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in  turn  enabled  integration  over  more  wave  cycles  before  the  computation  was  terminated 
due  to  suspected  wave  reflection  from  the  downstream  boundary23  (currently  up  to  G  wave 
cycles);  more  points  in  the  radial  direction  (currently  approximate!}'  100)  and  a  better 
distribution  of  mesh  points  next  to  the  wall  (the  first  mesh  point  a  distance  of  only  4.2  /mi 
away  from  the  wall);  and  finally,  to  help  overcome  the  first  order  accuracy  of  the  MINT 
code,  twice  as  many  time  steps  per  wave  cycle  are  utilized  (400  time  steps  compared  to 
200  previously).  In  addition.  Scientific  Research  Associates  of  Hartford  C'T,  graciously 
made  available  to  the  authors  a  partially  vectorized  version  of  the  MINT  code  which  is 
both  faster  (approximately  3.S  times  faster)  than  the  scalar  code  and  due  to  improved 
code  architecture,  requires  less  storage  for  the  identical  grid  size.  The  computational  mesh 
for  the  flow  refraction  study  utilizes  250  mesh  points  uniformly  distributed  in  the  axial 
direction  and  up  to  100  points  in  the  radial  direction.  Very  fine  mesh  distribution  is  utilized 
near  the  wall,  with  an  exponential  growth  of  mesh  spacing  toward  the  centerline. 

NUMERICAL  RESULTS 

The  first  test  case  utilized  a  frequency  of  6000  Hz  and  centerline  Mach  number  of  0.1. 
The  time  evolution  of  the  radial  distribution  of  acoustic  pressure  is  shown  in  Fig.  2a.  The 
initial  plane  wave  is  rapidly  distorted  as  the  wave  propagates  downstream.  After  600  time 
steps,  the  amplitude  of  the  acoustic  pressure  at  the  wall  is  39.5%  larger  than  at  centerline. 
The  corresponding  time  evolution  of  the  radial  distribution  of  the  axial  acoustic  velocity 
is  shown  in  Fig.  2b.  The  initial  plane  wave  is  distorted  in  time.  After  GOO  time  steps, 
the  amplitude  of  axial  velocity  oscillations  near  the  wall  is  31%  larger  than  at  centerline. 
It  should  be  pointed  out  that  the  amplitudes  of  both  the  acoustic  pressure  and  the  axial 
acoustic  velocity  were  decreased  as  the  wave  propagated  downstream. 

These  calculations  were  performed  both  on  the  Cray  1-S  computer  and  the  Cray  X- 
MP  computer.  The  calculations  on  the  Cray  X-MP  computer  were  performed  with  a 
significant  increase  of  mesh  points  per  wave  cycle  and  mesh  points  in  the  radial  direction, 
combined  with  a  significant  decrease  in  the  integration  time  step.  Nevertheless,  the  repeat 
calculations  resulted  in  pressure  and  axial  acoustic  velocity  value's  (and  thus,  based  on 
t ho  definition  of  acoustic  energy  as  p\v;  acoustic  energy)  that  were  not  significantly  differ¬ 
ent  (within  a  fe  v  percent)  from  the  results  obtained  by  utilizing  the  Cray  1-S  computer. 
Phis  observation  does  not  imply  that  artificial  nunfuical  dissipation  effects  had  been  com¬ 
pletely  eliminated.  Nevertheless,  it  suggests  that  nunn  rical  dissipation  constitutes  only  an 
insignificant  part  of  the  total  , acoustic  energy  dissipated  as  the  ware  traverses  the  length 
of  the  computational  domain. 
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Several  other  tests  with  different  frequencies  and  centerline  Mach  numbers  indicated 
that  acoustic  refraction  effects:  (a)  increase  with  frequency;  (b)  increase  with  centerline 
Mach  number;  and  (c)  are  higher  for  upstream  wave  propagation  than  for  downstream 
wave  propagation.  In  this  connection  it  should  be  mentioned  that  linear  theory  results 
indicated  the  same  trends.  Nevertheless,  the  large  refraction  effects  predicted  by  linear 
analysis  are.  as  expected,  significantly  higher  than  those  obtained  via  a  solution  of  the 
Navier-Stokes  equations. 

Based  on  the  classical  solution  for  acoustic  wave  propagation  in  the  presence  of  mean 
flow,  the  acoustic  boundary  layer  thickness  was  calculated  to  be  330  //in  for  a  frequency 
of  1000  Hz  and  132  /an  for  G000  Hz,  for  a  fully  developed  pipe  flow  with  centerline  Mach 
number  equal  to  0.1.  The  additional  (over  previous  calculations)  number  of  mesh  points 
distributed  in  the  radial  direction  enabled  better  resolution  of  the  acoustic  boundary  layer 
and  better  simulation  of  processes  that  describe  transfer  of  information  in  the  radial  direc¬ 
tion,  especially  inside  the  acoustic  boundary  layer.  A  wealth  of  information  was  obtained 
and  several  processes  have  been  observed.  Due  to  space  and  time  limitations,  only  a  few 
of  these  will  be  highlighted.  Discussion  of  other  results  and  further  analysis  of  the  results 
highlighted  here  will  be  deferred  to  future  publications. 

Figures  3  and  4  show  the  radial  distribution  of  the  acoustic  pressure  and  axial  velocity 
at  an  axial  station  located  1.32  m  downstream  from  the  left  boundary,  for  the  complete 
radius  of  the  tube,  and  for  an  narrow  zone  near  the  wall,  respectively.  This  test  was 
conducted  at  a  longitudinal  frequency  of  oscillations  of  1000  Hz.  The  radial  distance  is 
nondimensionalized  by  the  radius  (0  05  m),  with  the  origin  at  the  centerline.  The  results 
in  Figs.  3a  and  3b  show  typical  radial  distribution  of  acoustic  pressure  and  axial  acoustic 
velocity.  These  figures  show  the  radial  distribution  of  acoustic  pressure  and  axial  velocity 
at  time  equal  5.32  milliseconds  (abbreviated  henceforth  to  ms)  and  t=  5.72  ms,  when  the 
acoustic  wave  is  near  its  maximum  and  minimum,  respectively,  at  this  location.  Notice 
that  the  change  of  the  acoustic  pressure  across  the  tube  is  relatively  small  at  this  time 
in  the  cycle.  In  this  connection  it  should  be  mentioned  that  in  a  previous  study23  it  was 
shown  that  the  maximum  variation  of  acoustic  pressure  across  the  tube  cross-section  for  a 
downstream  propagating  wave  at  1000  Hz  was  less  than  1%.  Figures  4a  through  4f  describe 
the  transition  of  the  acoustic  pressure  and  axial  acoustic  velocity  from  the  positive  to  the 
negative  sign  of  waveform,  where  the  positive  waveform  corresponds  to  a  positive  acoustic 
pressure  and  axial  acoustic  velocity  (downstream  directed)  and  a  negative  wave  corresponds 
to  negative  acoustic  pressure  and  axial  acoustic  velocity  (upstream  directed).  At  t=4.75  ms 
(shown  in  Fig.  4a)  the  acoustic  wave  just  started  decaying  from  its  maximum  value. 


Notice  the  increase  of  acoustic  pressure  across  the  acoustic  boundary  layer  (approximately 
150  pm  at  this  time)  that  results  from  acoustic  refraction,  and  the  overshoot  of  acoustic 
axial  velocity  at  the  edge  of .  the  acoustic  boundary  layer  (Richardson’s  annular  effect). 
Figure  4b  shows  that  while  the  amplitudes  of  the  acoustic  pressure  and  axial  velocities 
are  decreasing,  a  further  decrease  is  observed  near  the  wall  (within  the  acoustic  boundary 
layer).  This  phenomenon  is  further  enhanced  at  t=4.96  ms  (Fig.  4c),  at  which  time  the 
axial  acoustic  velocity  near  the  wall  becomes  negative.  At  t=5.03  ms  (Fig.  4d)  the  acoustic 
pressure  across  the  tube  is  still  completely  positive,  while  the  axial  acoustic  velocity  is 
positive  everywhere,  except  near  the  wall.  At  this  time  the  thickness  of  the  newly  formed 
acoustic  boundary  layer  is  only  GS  pm.  It  should  be  noticed  that  simultaneously  with 
the  formation  of  a  new  acoustic  boundary  layer  near  the  wall  (with  a  negative  velocity 
sign),  the  previous  acoustic  boundary  layer  (with  a  positive  sign  and  its  own  overshoot) 
is  propagating  further  into  the  core  flow,  with  a  thickness  of  approximately  500  pm  (at 
t=5.02  ms).  The  acoustic  pressure  and  axial  acoustic  velocity  become  negative  only  at 
t=5.05  ms,  as  shown  in  Fig.  4e.  At  a.  later  time  (5.25  ms,  shown  in  Fig.  4f)  the  acoustic 
pressure  profile  across  the  acoustic  boundary  layer  and  across  the  tube  diameter  decrease 
in  value,  in  accordance  with  refraction  theory  for  upstream  wave  propagation. 

The  transition  from  a.  positive  waveform  to  a  negative  waveform,  shown  in  Figs.  4a 
through  4f,  demonstrates  that  acoustic  refraction  effects  result  in  increase  of  the  absolute 
value  of  the  acoustic  pressure  with  distance  from  the  centerline,  independent  of  the  sign 
of  the  wave.  When  the  wave  is  positive,  the  acoustic  pressure  value  near  the  wall  is  higher 
than  the  acoustic  pressure  value  at  centerline.  When  the  wave  is  negative  negative,  the 
acoustic  pressure  value  near  the  wall  is,  again,  higher  (in  absolute  value)  than  the  acoustic 
pressure  at  centerline.  Another  phenomenon  which  is  observed  is  the  time  variation  of  the 
acoustic  boundary  layer  thickness.  When  the  wave  is  near  its  maximum  value  (Fig.  4a), 
acoustic  boundary  layer  thickness  is  approximately  150  pm.  As  the  wave  changes  its  sign, 
the  acoustic  boundary  layer  thickness  starts  growing  from  a  zero  value  to  GS  pm  (Fig.  4d) 
and  finally  (Fig.  4f)  reaches  the  same  thickness  as  when  the  wave  was  at  its  maximum 
value  (i.e.,  about  150  pm).  Richardson’s  annular  effect  at  the  edge  of  the  acoustic  boundary 
layer,  which  is  somewhat  diminished  in  Figs.  3  due  to  scaling,  is  better  observed  in  Figs.  4. 
In  this  connection  it  should  be  mentioned  that  the  transition  from  a  negative  to  a  positive 
wave  is  of  similar  nature,  i.e.,  the  waveform  transition  occurs  for  the  first  time  near  the 
wall. 

Another  transition  of  interest  is  that  of  the  radial  acoustic  velocity  from  a  velocity 
directed  away  from  the  wall  (negative  sign)  to  a  velocity  directed  toward  the  wall  (positive 
sign).  Before  analyzing  these  results  it  should  be  mentioned  that  the  analysis  predicts 


that  the  radial  acoustic  velocity  leads  the  acoustic  pressure  and  axial  acoustic  velocity  by 
90  degrees  and  that  the  maximum  value  of  the  radial  acoustic  velocity  is  approximately 
two  orders  of  magnitude  smaller  than  the  maximum  value  of  the  axial  acoustic  velocity. 
The  transition  of  the  radial  acoustic  velocity  is  demonstrated  utilizing  results  obtained 
for  an  acoustic  wave  propagation  downstream  (centerline  Mach  number  equal  to  0.1)  at 
a  frequency  of  G000  Hz.  Figures  5a  through  5e  show  the  radial  distribution  of  the  radial 
acoustic  velocity  for  several  times  at  a  tube  cross  section  located  0.15S  m  downstream  from 
the  left  boundary.  At  t =0.5  ms  the  radial  acoustic  velocity  has  a  minimum  at  approxi¬ 
mately  0.G0  radii.  At  later  times  the  wave  is  distorted  with  the  minimum  approaching  the 
wall  while  the  positive  component  of  the  radial  acoustic  velocity  starts  growing  near  the 
centerline,  at  approximately  0.3  radii  (as  shown  in  Figs.  5b  through  5c).  During  this  time, 
the  acoustic  pressure  and  axial  acoustic  velocity  arc  approaching  their  minima  value,  which 
they  reach  at  t=0.5G  ms.  At  this  time  the  radial  acoustic  velocity  is  completely  positive 
with  an  almost  perfect  symmetric  distribution  across  the  tube  cross-section  (Fig.  5f).  The 
physical  implication  of  the  scenario  depicted  in  this  set  of  figures  is  that  during  the  half 
wave  cycle  at  which  the  acoustic  pressure  and  axial  acoustic  velocity  transition  from  a 
maximum  value  to  their  minimum  value,  the  flow  is  pushed  away  from  the  wall  towards 
the  centerline  (i.e.  negative  radial  acoustic  velocity).  During  the  other  half  wave  cycle, 
when  the  radial  acoustic  velocity  is  positive  and  the  acoustic  pressure  and  axial  acoustic 
velocity  transition  from  a  minima  values  to  their  maxima,  values,  the  flow  is  pushed  toward 
the  wall.  It  is  postulated  here  that  acoustic  refraction  phenomenon  is  strongly  influenced 
by  the  direction  of  the  radial  acoustic  velocity.  As  the  mass  of  the  fluid  is  driven  toward 
the  wall  by  the  positive  radial  acoustic  velocity,  the  acoustic  pressure  at  the  wall  becomes 
higher  than  the  acoustic  pressure  at  the  centerline.  When  the  flow  is  driven  toward  the 
centerline,  the  reverse  accurs.  Since  radial  velocity  values  are  approximately  two  orders 
of  magnitude  lower  than  axial  acoustic  velocity  values,  the  effect  of  refraction  is  small.  In 
this  connection  it  should  be  mentioned  that  for  downstream  wave  propagation  (for  either 
1000  Hz  or  G000  Hz),  it  was  found  that  the  maximum  value  of  the  radial  acoustic  velocity 
is  approximately  twice  as  large  as  the  minimum  value. 

To  better  correlate  the  results  obtained  with  results  available  in  the  literature18,  results 
that  are  typically  in  the  form  of  rms  values,  a  computation  of  the  acoustic  values  obtained 
at  1000  Hz  were  performed.  A  typical  result  is  shown  at  an  axial  station  located  0.515  m 
downstream  of  the  left  boundary.  Figure  6  shows  the  near  wall  radial  distribution  of  the 
acoustic  pressure  and  the  axial  and  radial  components  of  the  acoustic  velocity.  The  axial 
acoustic  velocity  figure  demonstrates  Richardson's  annular  effect .  The  rms  calculated  value 
for  the  acoustic  boundary  layer  thickness  is  approximately  ISO  //in.  in  good  agreement 


with  approximate  analysis  predictions20  and  experimental  trends18.  In  this  connection 
it  should  be  mentioned  that,  to  the  best  of  the  authors  knowledge  these  are  the  first 
solutions  demonstrating  cither  Richardson’s  annular  effect  or  refraction  effects  that  are 
obtained  through  solution  of  the  complete  Navier-Stokes  equations. 

As  noted  in  the  introduction,  it  was  expected  that  the  part  of  the  acoustic  energy  of 
the  waveform  will  be  converted  to  mean  flow  energy.  It  was  also  expected  that  the  planar 
wave  introduced  at  the  left  boundary  will  exhibit  an  increasingly  larger  mean  pressure  and 
axial  velocity  shifts  (indicating  increase  of  the  mean  pressure  and  axial  velocity  values, 
respectively)  as  the  wave  propagated  downstream.  Only  small  mean  pressure  or  axial 
velocity  shifts  is  observed  close  to  the  left  boundary.  The  phenomenon  becomes  more 
significant  as  the  wave  propagates  downstream.  Figures  7a  through  7c  show  the  time 
evolution  of  the  acoustic  pressure,  axial  acoustic  velocity  and  acoustic  energy  for  axial 
stations  35,  95,  and  155  located  on  the  centerline  at  distances  of  0.291,  0.806  and  1.321  m, 
respectively,  from  the  left  boundary.  It  is  shown  that  at  station  35  (Fig.  7a)  there  is 
virtually  no  mean  acoustic  pressure  or  velocity  shift  and  that  the  velocity  and  pressure 
are  in  phase.  The  situation  changes  as  the  wave  propagates  downstream.  At  station  95 
(Fig.  7b)  there  is  an  observed  mean  pressure  shift  of  approximtely  0.15%  and  a  mean 
velocity  shift  of  0.1%.  In  addition,  due  to  the  mean  pressure  and  axial  velocity  shifts, 
their  minima  are  slightly  out  of  phase,  as  observed  in  the  acoustic  energy  information 
(i.e.,  negative  acoustic  energy).  It  should  be  noted  that  since  the  definition  of  acoustic 
energy  was  adopted  from  plane  wave  theory,  negative  acoustic  energy  simply  implies  that 
the  acoustic  pressure  and  axial  acoustic  velocity  have  opposite  signs  at  a  given  instant  of 
time,  and  cannot  be  interpreted  as  acoustic  energy  production  or  dissipation.  At  axial 
station  155  (Fig.  7c)  the  mean  deviations  are  even  larger.  The  mean  pressure  is  shifted  by 
approximately  0.25%  and  the  mean  axial  velocity  by  0.2%.  The  phase  shift  between  these 
quantities  is  larger  than  before,  as  shown  in  the  acoustic  energy  plot.  While  these  numbers 
seem  small  they  correspond  to  a  significant  percentage  of  the  initial  acoustic  energy  in  the 
system.  A  disturbance  with  an  initial  pressure  amplitude  of  4.0%)  contains  only  4.92%  of 
the  total  energy.  Thus,  a  mean  pressure  shift  of  0.25%  and  a  mean  velocity  shift  of  0.2% 
represent  a  transfer  of  9.1%  of  the  initial  acoustic  energy  to  the  mean  flow  energy.  The 
increase  in  mean  pressure  remains  approximately  constant  across  the  tube  cross  section 
at  a  given  axial  location.  The  mean  axial  velocity  shift  increases  only  slightly  across 
most  of  the  tube  cross-section  and  starts  changing  considerably  only  within  approximately 
120  ftm  away  from  the  wall,  inside  the  acoustic  boundary  layer.  Within  this  zone  the  axial 
velocity  shift  increases  significantly  and  so  does,  consequently,  the  phase  shift  between  the 
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axial  acoustic  velocity  and  the  acoustic  pressure.  At  a  distance  of  60  fim  from  the  wall 
the  mean  velocity  shift  increase  to  0.32%  and  at  a  distance  of  17  pm  from  the  wall,  the 
mean  velocity  is  shifted  by  approximately  0.57%.  It  is  of  interest  to  notice  that  inside 
the  acoustic  boundary  layer  the  phase  differences  between  the  acoustic  pressure  and  axial 
acoustic  velocity  increase  toward  the  wall.  At  the  wall,  the  axial  acoustic  velocity  leads  the 
acoustic  pressure  by  approximately  30-40  degrees.  The  time  evolution  of  acoustic  pressure, 
axial  acoustic  velocity  and  acoustic  energy  at  a  point  located  4.72  /an  away  from  the  wall 
is  shown  in  Fig.  7d.  While  the  time  evolution  of  the  acoustic  pressure  is  almost  identical  to 
that  observed  at  the  centerline  (Fig.  7c),  with  an  identical  mean  pressure  shift  (0.25%),  the 
axial  acoustic  velocity  exhibits  a  nonlinear  waveform  with  a  mean  shift  of  approximately 
0.6%.  These  values  imply  that  approximately  1S%  of  the  initial  acoustic  energy  in  the 
system  had  been  converted  to  mean  flow  energy.  The  phase  shift  between  the  acoustic 
pressure  and  axial  acoustic  velocity  is  large  enough  to  obtain  negative  acoustic  energy 
values  (i.e.  the  acoustic  pressure  and  axial  velocity  have  different  signs)  over  a  significant 
number  of  time  steps. 

ACOUSTIC  WAVE  PROPAGATION 
WITH  NO  MEAN  FLOW 

The  computational  mesh  in  the  radial  direction  utilized  for  this  test  case  was  identical 
to  the  one  used  for  the  refraction  test  case.  Since  the  objective  was  to  integrate  the  solution 
for  only  three  wave  cycles  (as  opposed  to  5  or  G  waveevdes  in  the  refraction  study),  only 
132  mesh  points  were  used  in  the  axial  direction.  The  number  of  mesh  points  per  wave 
cycle  was  held  constant  at  44.  The  longitudinal  acoustic  frequency  was  1000  Hz. 

NUMERICAL  RESULTS 

After  propagating  downstream  for  approximately  three  wavelengths,  the  amplitude  of 
the  acoustic  pressure  was  reduced  by  about  12%,  while  the  amplitude  of  the  axial  acoustic 
velocity  was  reduced  by  about  10%.  This  corresponds  to  roughly  8%  acoustic  energy 
reduction  over  three  wavelengths. 

As  pointed  out  in  the  introduction,  previous  experimental15'18  and  theoretical21  inves¬ 
tigations  demonstrated  that  when  an  acoustic  wave  propagates  in  a  hard  wall  tube,  with 
or  without  mean  sheared  flow,  there  is  an  excitation  of  a  mean  flow  component  at  the  edge 
of  the  acoustic  boundary  layer  (i.e.  Richardson  s  annular  effect).  The  energy  contained  in 
the  newly  developed  mean  flow  has  to  lie  taken  out  of  the  oscillatory  field.  This  transfer 
of  energy  from  the  acoustic  to  the  mean  field  by  excitation  of  mean  flow  field  velocity, 
often  referred  to  as  acoustic  streaming,  results  in  net  energy  gain  to  the  mean  flow  in  the 
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chamber  (that  may  or  may  not  exist  there*  initially)  and  a  net  loss  to  the  oscillatory  field 
in  the  chamber. 
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Figures  $a  through  Sd  show  the*  radial  distribution  of  tin*  acoustic  pressure,  axial  and 
radial  acoustic  velocities  across  the  tube  cross-section  and  expanded  views  near  the  wall 
at  an  axial  station  located  a  distance  of  0.37S  m  downstream  from  the  left  boundary. 
Figures  Sa  and  Sb  show  the  radial  variation  at  the  time  when  the  acoustic  wave  is  near  its 
minimum  at  this  location  (t  — l.SS  ms),  while  Figs. Sc  and  Sd  show  the  radial  distribution 
half  a  wave  cycle  later,  at  t =2.3 1  ms.  At  this  time  the  acoustic  wave  is  near  its  maximum 
at  this  location.  The  radial  distribution  of  acoustic  pressure  demonstrates  small  refraction 
('fleet.  As  mentioned  previously,  when  sound  waves  propagate  in  a  coexisting  sheared 
mean  How.  the  effect  of  velocity  gradients  is  to  refract  the  sound.  Since  the  solution  was 
initiated  with  no  mean  flow,  the  small  refraction  effect  observed  indicates  excitation  of  a 
sheared  mean  flow  by  tin*  acoustic  field.  The  steeper  increase  (in  absolute  value)  of  acoustic 
pressure  inside  the  acoustic  boundary  layer  is  attributed  to  the  larger  axial  acoustic  velocity 
gradient  at  this  location. 

The  acoustic  boundary  layer  thickness  is  observed  to  vary  in  time  during  the  cycle 
with  a  minimum  value  at  the  time  the  acoustic  wave  changes  its  sign  and  a  maximum  value 
at  the  time  the  acoustic  wave  is  near  its  maximum  or  minimum  value.  Figures  9a  through 
9d  show  the  radial  variation  of  the  axial  acoustic  velocity  near  the  wall  at  four  time  instants 
during  the  transition  of  the  wave  from  a  negative  waveform  to  a  positive  waveform.  Figure 
9a  shows  the  radial  distribution  of  the  axial  acoustic  velocity  at  t  =  2 . 0 1  ms,  just  before 
the  wave  transitions  from  a  negative  to  a  positive  waveform.  The*  acoustic  boundary  layer 
thickness  is  203  /an.  Figure  9b  (at  t=  2.07  ms)  shows  that  a  new  downstrein  directed 
acoustic  boundary  layer  has  been  formed  near  the  wall,  while  the  axial  acoustic  velocity 
everywhere  else  is  still  directed  upstream  (negative  sign).  Notice  that  the  overshoot  at  the 
edge  of  the  acoustic  boundary  layer  is  propagating  toward  the  core  flow  as  the  acoustic 
boundary  layer  is  growing.  At  t=2.09  ms  (Fig.  9c)  the  downstream  directed  acoustic 
boundary  layer  continues  to  grow  in  value  and  expands  deeper  into  the  core  flow,  while 
the  edge  of  the  upstream  directed  boundary  layer  is  now  located  320  /nil  from  the  wall. 
Finally,  at  t  =  2 . 1 2  ms  (Fig.  9d)  the  axial  acoustic  velocity  across  the  complete  tube  is 
positive,  with  the  large  overshoot  at  the  edge  of  the  acoustic  boundary  layer,  located  a 
distance  of  03  //m  from  the  wall. 

The  mis  calculated  values  at  this  axial  location  are  shown  in  Fig.  10.  The  figure  shows 
the  acoustic  boundary  layer  to  have  an  average  thicknesss  of  230  /tin.  while  Richardson22 
predicts  a  value  of  330  //in.  Since  it  is  suspected  that  Richardson  could  not  account  for  the 
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type  of  variation  in  acoustic  boundary  layer  thickness  depicted  in  Figs.  9a  through  9d,  it 
is  postulated  that  the  value  of  330  pm  is  the  maximum  thickness  of  the- acoustic  boundary 
layer.  This  maximum  value  will  yield  an  rms  value  of  233  //m.  in  good  agreement  with 
the  computational  predictions.  In  this  connection  it  should  be  mentioned  that  the  rms 
acoustic  boundary  layer  thickness  with  no  flow  is  thicker  than  acoustic  boundary  layer 
thickness  with  mean  flow  (230  pm  vs  ISO  /mi),  in  agreement  with  theoretical  predictions21 
and  experimental  data15  that  show  boundary  layer  thickness  to  decrease  with  increase  of 
Reynolds  number. 

FLOW  TURNING  STUDY 

The  only  test  case  conducted  under  this  study  was  acoustic  wave  propagation  down¬ 
stream  in  a  tube  where  the  mean  flow  is  injected  into  the  tube  through  its  lateral  boundary 
(i.e.  tube  walls).  The  driven  longitudinal  acoustic  frequency  of  oscillations  was  1000  Hz. 
The  geometry  used  for  this  test  case  is  shown  in  Fig.  11.  The  length  of  the  injection  zone 
is  1.167  m,  a  distance  that  corresponds  to  approximately  three  wavelengths  at  a  frequency 
of  1000  Hz.  There  are  137  points  distributed  uniformly  in  the  axial  direction  along  the 
blowing  zone;  yielding  approximately  45  points  per  wave  length.  The  injection  Mach  num¬ 
ber  is  equal  to  0.0021.  After  a  steady  state  flow  solution  was  obtained,  the  pressure  and 
axial  velocity  at  the  upstream  (left)  boundary  were  perturbed  sinusoidally  to  produce  a 
system  of  downstream  traveling  waves  in  the  computational  domain. 

NUMERICAL  RESULTS 

Figures  12a,  12b.  13a,  and  13b  show  the  radial  distribution  of  acoustic  pressure,  axial 
acoustic  velocity  and  radial  acoustic  velocity  at  two  instants  during  a  single  wave  cycle. 
Data  are  shown  in  the  two  figures  for  the  complete*  radius  and  for  an  expanded  view  of 
the  radial  variation  near  the  blowing  wall.  These  results  were  obtained  at  an  axial  station 
located  1.004  m  downstream  from  the  left  boundary  and  0.1C3  m  upstream  from  the  end  of 
the  blowing  zone.  Figures  12a  and  13a  show  the  radial  variation  of  the  acoustic  properties 
at  t=4.24  ms.  At  this  time  the  acoustic  pressure  and  the  axial  acoustic  velocity  are  about 
to  reach  their  minima  values  while  the  radial  acoustic  velocity  is  about  to  transition  from 
a  velocity  directed  toward  the  centerline  (negative  sign)  to  a  velocity  directed  toward  the 
wall  (positive  sign).  It  is  noted  that  for  acoustic  wave  propagation  downstream  with  wall 
blowing,  the  radial  acoustic  velocity  leads  the  acoustic  pressure  by  90  degrees,  as  was  ob¬ 
served  for  the  previous  two  studies.  The  figures  demonstrate  that  due  to  blowing  from  the 
wall,  the  acoustic  boundary  layer  is  significantly  thicker  than  for  acoustic  wave  propagation 
with  or  without  mean  flow.  The  acoustic  boundary  layer  thickness  varies  between  260  to 
approximately  500  //in,  as  compared  to  rms  values  of  ISO  //m  for  the  refraction  test  case 


15 


1 


1 


I 

I 


and  250  /nn  for  the  no  mean  flow  case.  The  axial  acoustic  velocity  overshoot  at  the  edge 
of  the  acoustic  boundary  layer  (Richardson’s  effect  )  is  observed  again.  However,  for  the 
wall  blowing  test  case  this  phenomenon  is  overshadowed  by  formation  of  an  adjacent  large 
positive  (i.e.,  downstream  directed)  increase  in  the  axial  acoustic  velocity.  This  positive 
increase  has  its  maximum  value  at  a  distance  of  approximately  7G0  /mi  away  from  the  in¬ 
jection  wall.  The  axial  acoustic  velocity  decays  from  this  maximum  overshoot  value  to  its 
centerline  value  only  at  a  distance  of  approximately  5000-G000  // m  away  from  blowing  sur¬ 
face.  The  value  of  the  axial  velocity  at  the  maximum  is  0.115  nondimensional  units  larger 
than  the  centerline  value.  In  this  connection  it  should  be  recalled  that  acoustic  properties 
are  defined  here  as  the  difference  between  the  time  dependent  and  the  steady  properties. 
Thus,  they  include,  in  addition  to  true  acoustic,  properties,  the  temporal  changes  of  the 
mean  flow  field. 

Figures  12b  and  131)  show  the  radial  distribution  of  acoustic  pressure,  axial,  and  radial 
acoustic  velocities  at  t=4.G7  ms;  a  difference  of  approximately  half  a  wave  cycle  from  the 
time  for  which  the  results  presented  in  Figs.  12a  and  13a  were  obtained.  At  this  time 
the  acoustic  pressure  and  acoustic  axial  velocity  are  about  to  reach  their  maxima  values 
while  the  radial  acoustic  velocity  is  about  to  change  from  a  wall  directed  to  a  centerline 
directed  (i.e.  from  a  positive  to  a  negative  value).  The  variation  of  acoustic  pressure 
across  the  acoustic  boundary  layer  (due  to  refraction)  is  demonstrated  in  Fig.  13b.  The 
overshoot  of  the  axial  acoustic  velocity  at  the  edge  of  the  acoustic  boundary  layer  is  not 
immediately  obvious  in  this  figure,  since  the  overshoot  of  the  axial  acoustic  velocity  and  the 
larger  value  axial  acoustic  velocity  positive  overshoot  that  was  also  described  previously, 
are  now  merged  (since  both  are  positive.)  It  is  important  to  notice  that  the  value  of  the 
axial  acoustic  velocity  at  the  maximum  of  the  "large  positive  overshoot”  is  again  0.115 
nondimensional  units  higher  than  the  value  of  the  axial  acoustic  velocity  at  the  centerline. 
An  examination  of  several  other  radial  distributions  of  the  axial  acoustic  velocity  at  other 
times  reveals  that  this  large  overshoot  exists  at  any  instant.  Moreover,  the  maximum 
axial  acoustic  velocity  value  is  consistently  0.115  nondimensional  units  higher  than  the 
value  of  the  centerline  axial  acoustic  velocity.  In  addition,  although  the  centerline  axial 
velocity  changes  signs  every  half  wave  cycle, this  large  scab'  axial  velocity  increase  is  always 
oriented  downstream  (i.e..  in  the  plus  directon).  This  indicates  an  addition  to  the  mean 
fixial  velocity,  and  thus,  an  addition  to  the  mean  flow  energy,  in  a  zone  extending  from  the 
edge  of  the  acoustic  boundary  layer  to  a  distance'  of  approximately  5000-G000  //in  away  from 
the  injection  wall.  The  value  and  the  location  of  the  maximum  value  of  this  mean  velocity 
shift  (referred  to  henceforth  as  axial  velcoity  mean,  or  dc,  shift  )  varies  with  distance  from 
the  left  bound/try,  as  will  be  shown  later.  Nevertheless,  it  should  be  emphasized  that  at  a 
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given  location  in  the  tube,  the  value  of  the  axial  velcoity  dc  shift  is  constant,  independent 
of  the  value  or  the  direction  of  the  axial  acoustic  velocity  at  the  centerline.  To  the  best  of 
the  authors  knowledge,  this  phenomenon  has  never  been  reported  before. 

The  transition  from  a  positive  wave  to  a  negative  wave  that  was  examined  for  the 
wave  refraction  study  is  significantly  more  complex  for  this  case  and  involves  information 
exchange  between  the  radial  acoustic  velocity,  pressure  and  axial  acoustic  velocity.  A 
complete  description  of  this  transition  would  occupy  several  pages  of  discussion  and  figures 
and,  due  to  lack  of  space,  will  be  left  for  future  papers.  In  this  paper  only  the  transition  of 
axial  acoustic  velocity  near  the  surface  will  be  examined.  Figures  14a  through  14g  show  the 
radial  distribution  of  the  axial  acoustic  velocity  near  the  blowing  wall  at  several  instants 
of  time,  at  an  axial  station  located  1.004  m  from  the  left  boundary.  At  t=4.19  ms  (shown 
in  Fig.  14a),  when  the  axial  acoustic  velocity  is  near  its  minimum  value,  the  thickness 
of  the  acoustic  boundary  layer  is  approximately  250  /an,  and  the  overshoot  of  the  axial 
acoustic  velocity  at  the  edge  of  the  acoustic  boundary  layer  is  large.  As  the  acoustic 
pressure  and  axial  acoustic  velocity  values  increase  from  their  minima,  transitioning  to 
their  maxima  values,  the  overshoot  is  slightly  diminished  (as  shown  in  Fig  14b).  Figure 
14c  (at  t=4.40  ms)  shows  that  as  the  value  of  the  centerline  axial  acoustic  velocity  is 
further  increased  (though  still  negative),  the  axial  acoustic  velocity  near  the  wall  becomes 
positive.  An  intriguing  situation  is  further  developed  as  a  new  acoustic  boundary  layer  is 
developed  inside  the  old  acoustic  boundary  layer;  the  axial  acoustic  velocity  at  a  zone  near’ 
the  wall  is  directed  downstream  while  further  away  from  the  wall,  inside  the  old  acoustic 
boundary  layer  (with  a  thickness  of  410  pm),  as  well  as  in  the  core  flow,  the  velocity  is 
still  directed  upstream  (at  time  equal  4.43  ms  and  4.45  ms,  as  shown  in  Figs  14d  and  14e, 
respectively).  The  axial  acoustic  velocity  is  directed  downstream  across  the  complete  tube 
only  at  t=4.51  ms.  The  downstream  directed  acoustic  boundary  layer  thickness  now  starts 
growing  with  time.  The  overshoot  at  the  edge  of  the  acoustic  boundary  layer  reaches  a 
limiting  value  while  the  axial  acoustic  velocity  everywhere  else  increases.  Finally,  as  shown 
in  Figs.  14f  and  14g  (for  t=4.53  and  t=4.61  ms,  respectively),  the  large  axial  velocity  dc 
shift  merges  with  the  overshoot  at  the  edge  of  the  acoustic  boundary  layer. 

As  stated  earlier,  this  discussion  is  merely  an  indication  of  the  complexity  of  the 
processes  by  which  the  acoustic  signal  changes  its  sign.  For  a  full  understanding  of  this 
process,  the  acoustic  pressure  and  the  axial  and  radial  acoustic  velocities  need  to  be  ana¬ 
lyzed  simultaneously  at  a  region  very  near  the  wall  as  well  as  across  the  tube  cross  section. 
Nevertheless,  even  with  the  limited  amount  of  information  shown,  it  was  clearly  demon¬ 
strated  that,  the  transition  is  initiated  at  the  blowing  surface.  The  core  flow  transition 
occurs  only  after  the  transition  inside  the  acoustic  boundary  layer  had  been  completed. 


ROOT  MEAN  SQUARE  AND  TIME  AVERAGED  CALCULATIONS 

The  root  mean  square  calculations  were  expected  to  show  the  acoustic  boundary  layer 
thickness  and  Richardson  s  annular  effect.  The  radial  distribution  of  the  root  mean  square 
(rms)  calculated  values  of  the  acoustic  pressure  and  axial  and  radial  acoustic  velocities  are 
shown  in  Figs.  15a  and  lob.  for  the  complete  radius  and  for  an  enlarged  view  of  a  zone  near 
the  wall,  respectively,  at  a  station  located  0.377  in  from  the  left  boundary.  The  acoustic 
pressure  data  shows  a  small  positive  refraction  effect  across  the  tube  diameter,  a  small 
undershoot  at  the  edge  of  the  acoustic  boundary  layer  and  an  added  acoustic  pressure 
increase  across  the  acoustic  boundary  layer.  The  rms  axial  acoustic  velocity  information 
demonstrates  Richardson’s  annular  effect  with  a  significant  overshoot  at  the  edge  of  the 
acoustic  boundary  layer.  The  calculated  rms  thickness  of  the  acoustic  boundary  layer  is 
330  /mi.  It  is  of  interest  to  note  that  the  radial  acoustic  velocity  also  has  an  increase  at 
the  edge  of  the  acoustic  boundary  layer,  unlike  the  rms  calculated  value  for  the  mean  flow 
claculations  and  similar  to  the  results  obtained  for  no  mean  flow.  In  this  connection  it 
should  be  mentioned  that  the  radial  acoustic  velocity  values  for  the  wall  blowing  test  case 
are  significantly  larger  (at  maximum  or  minimum)  than  the  corresponding  values  obtained 
for  the  refraction  test  case,  due  to  mean  flow  injection  in  the  radial  direction.  In  addition, 
it  is  noted  that  the  maxima  values  of  the  radial  acoustic  wave  were  approximately  twice 
as  large  as  the  minima  values,  as  was  observed  for  the  flow  with  no  wall  injection. 

Since  mass  is  injected  into  the  computational  domain  only  through  the  forward  part 
of  the  tube,  downstream  of  the  blowing  zone  the  flow  reverts  slowly  to  the  same  flow 
pattern  observed  for  acoustic  wave  propagation  in  the  presence  of  mean  flow  but  no  wall 
blowing  (i.e.  the  refraction  test  case).  Figure  15c  shows  the  radial  distribution  of  acoustic 
pressure  and  axial  and  radial  acoustic  velocities  near  the  hard  wall  at  a  station  located 
0.112  m  downstream  from  the  end  of  the  blowing  zone  (1.279  m  from  the  left  boundary). 
The  acoustic  boundary  layer  shrank  to  approximately  half  the  thickness  predicted  for  the 
blowing  zone  (rms  vine  of  190  // m  vs.  330  /mi  in  the  blowing  zone),  and  the  undershoots  and 
overshoots  of  acoustic  pressure  and  radial  acoustic  velocity,  respectively,  at  the  edge  of  the 
acoustic  boundary  layer  that  were  observed  near  the  blowing  zone,  had  now  disappeared. 

A  potentially  significant  new  phenomenon  that  was  observed  during  the  analysis  of  the 
transient  results,  namely,  the  axial  velocity  dc  shift  that  started  at  the  edge  of  the  acoustic 
boundary  layer,  was  not  evident  in  the  rms  calculations.  It  is  a  dc  phenomenon  and  as 
such  cannot  be  evident  in  rms  calculations.  However,  this  type  of  phenomenon  should  be 
evident  in  time  averaged  calculations.  Figures  IGa  through  lGf.  show  the  radial  distribution 
of  the  axial  acoustic  velocity  at  six  axial  locations;  stations  31.  51.  71.  91,  111  and  131 
are  located  at  distances  of  0.257  m,  0.429  m.  0.G00  m,  0.772  m.  0.944  m  and  1.115  in 


from  the  left  boundary*  respectively.  All  stations  are  located  within  the  blowing  zone. 
The  mean  axial  velocity  shift  maximum  is  initially  located  at  a  distance  of  approximately 
2250  /mi  away  from  the  blowing  wall  (Fig.  IGa).  As  the  acoustic  wave  (initiated  at  the  left 
boundary)  propagates  downstream  over  the  blowing  zone,  the  location  of  the  time  averaged 
mean  velocity  shift  maximum  is  approaching  the  wall.  At  axial  location  131  (Fig.  lGf), 
the  maximum  is  located  approximately  GGO  /an  away  from  the  wall.  In  addition,  as  the 
acoustic  wave  propagates  downstream  along  the  blowing  zone,  the  maximum  value  of  the 
axial  velocity  clc  shift  increases.  At  station  31  this  maximum  value  is  0.00G  nonclimensional 
units  higher  than  the  centerline  time  averaged  acoustic  axial  velocity  (Fig.  IGa).  0.015  at 
station  51  (Fig.  1Gb)  0.026  at  station  71  (Fig.  lGc),  0.035  at  station  91  (Fig.  16d),  0.042 
at  station  111  (Fig.  lGe)  and  0.04S  at  station  131  (Fig.  16f).  In  addition  to  this  important 
increase  in  the  maximum  value  of  the  axial  velocity  dc  shift,  these  figures  also  demonstrate 
another  very  important  fact;  namely,  the  time  averaged  axial  acoustic  velocity  at  most 
points  has  a  positive  value.  This  implies  that  acoustic  energy  had  been  converted  to  mean 
flow  energy  even  at  the  centerline;  again  demonstrating  acoustic  streaming  effects. 

The  results  obtained  for  the  wave  refraction  study  demonstrated  mean  pressure  and 
axial  velocity  shifts,  indicating  transfer  of  energy  from  the  acoustic  field  to  the  mean  flow 
field.  A  similar  investigation  was  conducted  for  the  acoustic  wave  propagation  in  the 
presence  of  wall  injection  by  analyzing  both  the  spatial  and  temporal  evolution  of  the 
acoustic  waveform  in  the  chamber.  Figures  17a  through  17c  show  the  spatial  evolution  of 
acoustic  pressure,  axial  acoustic  velocity  and  acoustic  energy,  at  t=4.20  ms,  as  a  function  of 
axial  distance  from  the  left  boundary,  for  radial  surfaces  located  on  the  centerline.  735  /an 
and  118  pm  away  from  the  wall,  respectively.  The  spatial  decay  of  both  acoustic  pressure 
and  axial  acoustic  velocity  as  well  as  the  growth  of  the  mean  pressure  and  axial  velocity 
shifts  are  shown  at  all  radial  surfaces.  It  should  be  noticed  that  the  disturbances  observed 
in  the  axial  acoustic  velocity  spatial  evolution  figures  for  radial  surfaces  near  the  wall  (Fig. 
17b  and  17c)  result  due  to  the  adjustment  of  the  flow  near  the  wall  from  a  flow  over  a 
blowing  wall  to  a  flow  over  a  hard  wall  (the  blowing  zone  ends  1.16  m  downstream  of  the 
left  boundary).  The  continuous  growth  of  the  axial  acoustic  velocity  (shown  in  Fig.  17b) 
is  another  demonstration  of  the  continuous  spatial  growth  of  the  axial  velocity  dc  shift 
discussed  in  the  previous  paragraph,  since  this  radial  surface  is  located  within  the  axial 
velocity  dc  shift  zone. 

Analysis  of  the  time  evolution  of  the  acoustic  pressure  and  axial  acoustic  velocity  at 
several  mesh  points  demonstrates  that:  (a)  the  mean  pressure  shift  increases  with  distance 
from  the  left  boundary  and  is  approximately  constant  across  the  tube  cross  section  (at 


any  given  axial  location);  (b)  the  mean  shift  of  axial  velocity  increases  both  in  the  axial 
and  radial  directions;  and  (c)  the  de  shifts  of  both  pressure  and  axial  acoustic  velocity 
result  in  the  growth  of  the  negative  component  of  the  acoustic  energy  (i.e..  the  number  of 
time  instants  in  which  the  acoustic  pressure  and  axial  velocity  have  opposite  signs)  both 
in  the  axial  and  radial  directions.  Figures  ISa  through  ISc  show  the  time  evolution  of 
axial  acoustic  velocity  tit  a  tube  cross-section  located  1.059  m  downstream  from  the  left 
boundary;  the  throe  stations  are  located  tit  the  centerline.  024  // m  from  the  blowing  surface 
(close  to  the  maximum  of  the  axial  velocity  de  shift)  and  90  //in  from  the  blowing  surface 
(within  the  acoustic  boundary  layer).  The  mean  pressure  shift  at  this  axial  location  is 
approximately  0.23%.  The  mean  axial  velocity  shift  is  only  0.1%  at  the  centerline  and 
grows  toward  the  wall.  Figure  ISb  shows  that  the  shift  in  axial  acoustic  velocity  still 
grows  after  one  wave  cycle.  Based  on  experience  gained  in  analyzing  results  obtained  at 
points  at  the  same  radius  but  closer  to  the  left  boundary,  points  that  experienced  4  or 
5  wave  cycles,  it  was  concluded  that  the  axial  velocity  mean  shift  stops  growing  after 
approximately  two  wave  cycles.  After  two  wave  cycles  the  mean  axial  velocity  shift  at  this 
location  is  approximately  0.9%.  These  values  of  mean  pressure  and  axial  velocity  shifts 
imply  that  at  this  location  approximately  23  percent  of  the  initial  acoustic  energy  had 
been  converted  to  mean  flow  energy.  Inside  the  boundary  layer  (Fig.  ISc)  the  mean  axial 
velocity  shift  is  approximately  0.5%,  while  the  waveform  is  nonlinear. 

ACOUSTIC  ENERGY  CALCULATIONS 

Figure  19  shows  the  comparison  of  normalized  time  averaged  acoustic  energy  as  a 
function  of  axial  distance  from  the  left  boundary  for  three  calculations:  namely:  acoustic 
wave  propagaton  with  no  mean  flow,  acoustic  wave  propagation  in  the  presence  of  mean 
flow  (acoustic  refraction  study)  and  acoustic  wave  propagation  in  a  tub/'  in  which  the  mean 
flow  is  injected  through  the  walls  of  the  tube  (flow  turning  study).  All  calculations  were 
performed  for  an  acoustic  wave  frequency  of  1000  Hz.  The  time  averaging  calculations 
were  performed  for  a  distance  of  only  one  wavelength  for  the  110-mean-flow  results  and  for 
three  wavelengths  (the  length  of  the  blowing  zone)  for  both  the  refraction  and  wall  blowing 
test  cases. 

The  figure  shows  that  as  the  wave  propagates  downstream,  reduction  of  acoustic 
energy  is  higher  with  mean  sheared  flow  than  with  no  mean  flow  and  even  higher  with 
wall  injection.  Acoustic  energy  is  reduced  in  1 1  if*  no  mean  flow  case  due  to  both  artificial 
numerical  dissipation  and  transfer  of  acoustic  energy  to  the  mean  flow,  as  was  shown 
during  the  discussion  of  the  results.  Since  the  exact  definition  of  acoustic  energy  for 


viscous,  rotational  flows  has  not  been  formulated  yet,  attempts  to  quantify  acoustic  energy 
transfer  to  the  mean  flow,  and.  by  subtracting  these  losses  from  the  total  acoustic  energy 
loss  evaluate  the  numerical  errors,  will  be  left  for  future  analysis.  Nevertheless,  it  is 
shown  that  losses  for  the  no-mean-flow  test  are  significantly  lower  than  losses  for  either 
the  refraction  or  wall  injection  test  cases.  Higher  acoustic  energy  losses  for  the  refraction 
test  case  may  indicate  that:  (a)  more  acoustic  energy  had  been  converted  to  mean  flow 
energy,  as  observed  in  the  analysis  of  the  time  dependent  results;  and  (b)  acoustic  energy 
dissipation  may  increase  due  to  the  mean  sheared  flow.  Finally,  the  figure  shows  that 
when  the  mean  flow  is  injected  through  the  wall,  acoustic  energy  losses  increase.  As  was 
noted  previously,  significantly  larger  mean  pressure  and  axial  velocity  shifts  were  observed 
throughout  the  computational  domain  for  the  flow  turning  test  case  than  for  the  refraction 
test  case.  .  In  addition,  an  even  larger  mean  flow  enhancement  was  observed  at  the  edge 
of  the  acoustic  boundary  layer.  Presently  it  is  not  immediately  apparent  how  to  quantify 
acoustic  energy  conversion  to  mean  flow  energy  in  a  viscous,  rotational  field.  Nevertheless, 
it  is  demonstrated  that  wall  injection  resulted  in  a  significantly  larger  acoustic  energy  loss, 
part  of  which  was  converted  to  mean  flow  energy  while  the  rest  was  dissipated. 

CONCLUSIONS 

This  paper  reports  initial  results  obtained  in  a  numerical  investigation  of  acoustic 
energy  exchange  between  the  mean  and  acoustic  flow  fields  through  acoustic  refraction, 
flow  turning  and  acoustic  streaming.  The  time-  dependent  compressible  Navier-Stokes 
equations  were  solved  utilizing  an  implicit,  non-iterative  Linearized  Block  Implicit  scheme. 
Three  test  cases  were  investigated:  acoustic  propagation  in  a  tube  with  a  coexisting  sheared 
mean  flow  (referred  to  as  a  refraction  study),  acoustic  propagation  in  a  tube  with  no  mean 
flow,  and  acoustic  propagation  in  a  tube  where  the  coexisting  mean  flow  is  injected  into 
the  tube  through  part  of  its  lateral  boundary  (referred  to  as  flow  turning  study). 

The  results  of  the  refraction  study  demonstrated  the  effect  of  acoustic  refraction. 
Sevc  ■ral  test  cases  were  repeated  on  a  more  powerful  computer  (Cray  X-MP)  with  reduced 
spatial  and  temporal  errors.  The  results  indicated  that  acoustic  refraction  effects;  (a)  in¬ 
crease  with  frequency:  (1))  increase  with  Mach  number;  (c)  are  higher  for  upstream  wave 
propagation  than  for  downstream  wave  propagation:  and  (d)  are  significantly  lower  than 
predicted  by  linear  theory.  Since  acoustic  refraction  effects  have  the  potential  to  enhance 
both  pand  whear  the  propellant  and  since  the  propellant  transient  burn  rate  is  strongly  de¬ 
pendent  on  these  flow  variables  (through  pressure  and  velocity  coupling  effects),  jt  appears 
that  acoustic  refraction  effects  should  be  included  in  solid  motor  stability  analvses. 


The  results  obtained  in  all  the  three  studies  demonstrated  Richardson’s  annular  ef¬ 
fect.  The  magnitude  of  the  overshoot  at  the  edge  of  the  acoustic  boundary  layer  and  the 
boundary  layer  thickness  changed  in  time  during  the  wave  cycle  and  changed  from  one  test 
case  to  another.  The  acoustic  boundary  layer  thickness  was  largest  for  the  wall  blowing 
case,  as  should  be  expected.  The  boundary  layer  thickness  was  larger  for  the  no-mean-flow 
case  than  for  the  refraction  test  case  since  acoustic  boundary  layer  thickness  decreases  with 
increase  in  Reynolds  number. 

All  three  cases  demonstrated  acoustic  streaming  effects.  The  largest  mean  pressure 
and  axial  velocity  shifts  were  obtained  for  the  flow  turning  case.  The  mean  pressure  shift 
increased  with  downstream  propagation  along  the  tube  and  was  approximately  constant 
at  any  axial  cross-section.  The  mean  axial  velocity  shift  also  increased  with  distance  from 
the  left  boundary.  In  addition,  the  axial  velocity  shift  increased  significantly  within  the 
acoustic  boundary  layer.  For  the  flow  turning  test  case,  an  even  larger  mean  axial  velocity 
dc  shift  was  observed  within  a  zone  adjacent  to  the  acoustic  boundary  layer,  extending 
a  distance  of  up  to  0.1  radii  away  from  the  wall.  The  distance  between  the  wall  and  the 
location  of  the  maximum  value  of  this  mean  axial  velocity  shift  was  reduced  with  increase 
of  dist  ance  from  the  left  boundary  while,  simultaneously,  the  maximum  value  increased  in 
value.  Based  on  initial  calculations  it  appears  likely  that  a  significant  amount  of  acoustic 
energy  is  converted  to  mean  flow  field  energy  within  this  zone. 

Complex  flow  phenomena  occur  as  the  acoustic  field  at  a  given  location  transitions 
from  positive  acoustic  pressure  and  axial  acoustic  velocity  to  negative  values,  or  from  a 
radial  acoustic  velocity  directed  toward  the  wall  (positive)  to  a  radial  velocity  directed 
toward  the  centerline  (negative  velocity).  It.  was  shown  that  the  radial  acoustic  velocity- 
leads  the  acoustic  pressure  and  axial  acoustic  velocity  by  90  degrees  and  that  near  the 
wall,  the  axial  acoustic  velocity  leads  the  acoustic  pressure  by  approximately  30  degrees. 

The  results  demonstrated  that  acoustic  energy  is  diminished  with  distance  from  the 
inflow  boundary.  For  acoustic  wave*  propagation  with  no  mean  flow,  this  acoustic  energy 
decrease  results  from  to  both  acoustic  energy  transfer  to  the  mean  flow  and  artificial 
energy  dissipation  due  to  numerical  errors.  Acoustic  energy  losses  increased  significantly 
when  sheared  mean  flow  was  introduced,  as  acoustic  dissipation  increased  due  to  wave 
propagation  in  a  viscous. rotational  flow  field.  Finally,  when  the  mean  flow  was  injected 
into  the  tube  through  the  lateral  boundary,  flow  turning  effects  resulted  in  a  significant 
increase  of  acoustic  energy  losses.  Though  the  present  research  verified  the  existence  of 
flow  turning  loss  phenomenon,  more  theoretical  and  experimental  efforts  art'  needed  to  help 
tpiantify  flow  turning  loss  anti  the  several  flow  phenomena  observed  in  this  investigation. 
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Fig.  5.  Radial  Distribution  of  Radial  Acoustic  Velocity,  f=6000  Hz,  Z=0.158  m.  (a)  Time=0.50  ms.  (b)  Time=  0.53  ms. 
(c)  Time=0.542  ms.  (d)  Time=0.546  ms.  (e)Time=  0.55  ms.  (f)  Time=0.56  ms. 
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Fig.  6.  Radial  Distribution  of  Root-Mean-Squared  Acoustic  Pressure,  Axial  Acoustic  Velocity  and  Radial  Acoustic  Velocity. 
f=1000  Hz,  Z=0.515  m. 
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Fig.  9.  Radial  Distribution  of  Acoustic  Velocity,  No  Mean  Flow,  f=1000  Hi,  Z=0.378  m.  (a)  Time=2.01  ms.  (b)  Time=2.07 
ms.  (c)  Time=2.09  ms.  (d)  Time=2.12  ms. 
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Fig.  10.  Radial  Distribution  of  Root  Mean  Squared  Acoustic  Pressure,  Axial  acoustic  Velocity  and  Radial  Acoustic  Velocity, 
No  Mean  Flow,  f=1000  Hz. 
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Fig.  11.  Computational  Grid  Used  for  Flow  Turning  Study. 
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Fig.  15.  Radial  Distribution  of  Root  Mean  Squared  Acoustic  Pressure,  Axial  acoustic  Velocity  and  Radial  Acoustic  Velocity 
Flow  Turning  Study,  (a)  Z=0.377  ra.  R=0.0-1.0.  (b)  Z=0.377  m.  R=0.945-1.0.  (c)  Z=1.279  m.  R=0.945-1.0. 
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Fig.  16.  Radial  Distribution  of  Time  Averaged  Axial  Acoustic  Velocity,  Flow  Turning  Study,  (a)  Z=0.257  m.  (b)  Z=0.429  m. 
(c)  Z=0.600  m.  (d)  Z=0.772  m.  (e)  Z=0.994  m.  (f)  Z=1.115  m. 
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Fig.  17.  Spatial  Evolution  of  Acoustic  Pressure,  Axial  Acoustic  Velocity  and  Acoustic  Energy,  Flow  Turning  Study,  Time=4.20 
ms,  (a)  Centerline,  (b)  Distance  to  Wall  =  735  microns,  (c)  Distance  to  Wa!l=118  microns. 
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Fig.  18.  Time  Evolution  of  Axial  Acoustic  Velocity,  Flow  Turning  Study,  Z=  1.059  m. 
microns,  (c)  Distance  to  Wall=94  microns. 
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Fig.  19.  Comparison  of  The  Axial  Distribution  of  Normalized  Acoustic  Energy  for  Refraction,  No  Mean  Flow  and  Flow  Turning 
Studies,  f=1000  Hz. 
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